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We present a theory of the four-terminal conductance for the multi-channel tunneling barrier, which 
is based on the self-consistent solution of Shrodinger, Poisson and continuity equations. We derive 
new results for the case of a barrier embedded in a long wire with and without disorder. We also 
recover known expressions for the conductance of the barrier placed into a ballistic constriction. 
Our approach avoids a problematic use of two chemical potentials in the same system. 
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I. INTRODUCTION 

Tunneling conductance through a potential barrier has 
been the subject of intense study since the original Lan- 
dauer paper [|] in 1957 (see Ref. for reviews). It was 
recognized that the two-terminal measurements of the 
tunneling conductance, g, may be described by the Lan- 
dauer formula 



e \ - , 



(1) 



where the ti are the eigenvalues of a barrier transmis- 
sion matrix. The term "two-terminal" implies that the 
voltage drop is measured between the same two contacts, 
sufficiently far from the barrier, which are used to pass a 
current. Actually Landauer's original expression ||l|], for 
the single channel case, had the different form 
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where r is the reflection coefficient of the barrier, obey- 
ing |rp -I- |tp — 1. Equation ^ is intuitively appealing 
since it leads to zero resistance for the fully transparent 
{\t\ = 1) barrier. Only after the work of Thouless ||] 
and Imry ||^ was it realized that Eq. (H) implicitly relies 
on the ability to measure the voltage drop right across 
the tunneling barrier. One should thus assume the exis- 
tence of another pair of contacts used solely as potential 
(as opposed to current) probes - a Jour-terminal mea- 
surement in modern terminology. The difference for the 



single channel case, g 



Trfi/e^, was interpreted 



as an effective resistance of the perfect wire, which con- 
nects the barrier with the current contacts. 

The two-terminal result, Eq. (^, may be derived from 
the Kubo formula for non-interacting electrons |^-^ . Its 
general form is unchanged by the Hartree screening of 
the barrier Q, which will, however, alter the transmis- 
sion coefficients. 

Thouless j^] confirmed Eq. (|^) for an electron gas, 
interacting on the Hartree level, by imposing charge 
neutrality in the wire on long length scales. In fact. 



the four-terminal expression, Eq. (g), is the result of 
a self-consistent calculation of the charge density in- 
duced near the barrier. Although this charge does 
not alter the divergence-less current, it does determine 
the voltage drop across the barrier region. To illus- 
trate this point let us consider a long single-channel 
wire with an embedded barrier. Suppose that the wire 
is placed in a time-dependent, spatially uniform ex- 
ternal electric field E{t) = E'^^* cos tot. This can be 
done, for example, by bending the wire into a ring, 
which encloses a time-dependent magnetic flux. Far 
from the barrier the electron current corresponds to 
pure acceleration, it is position independent and given 
by i{t) = e^n/(maj)i?'^^* sinwt, where n is the elec- 
tron density. Close to the barrier, however, the cur- 
rent is space dependent and charges accumulate on both 
sides of the barrier. These induce an additional electric 
field which contributes to the voltage drop, V , across 
the barrier. Calculations given below lead to the result 
V = E^'-"^ s\TL{Lut)2vF\r\'^ I iWio). Taking into account 
that in one dimension n = 2mwi?/(7r?i), where t^F is the 
Fermi velocity, one recovers Eq. (^ for the conductance, 
in agreement with Thouless's conclusion. 

While the single-channel four-terminal expression, 
Eq. (H), is generally agreed upon, the generalization to 
the multi-channel case has been a subject of controversy. 
P. W. Anderson et. al. [|j proposed to add the conduc- 
tances of channels in parallel. Langrcth and Abrahams 
1^ assumed that the chemical potential, say on the left of 
the barrier, is determined by the mean chemical poten- 
tial of all left and right moving electrons. The difference 
between mean left and mean right chemical potential was 
identified with the voltage drop. Calculating the conduc- 
tance required solving a large system of equations. Their 
approach assumes the existence of strong independent 
equilibration among the left- and right-moving particles, 
and the absence of any cross-talk between the left and 
right movers. As was first remarked by Biittiker, Imry, 
Landauer and Pinhas (BILP) ||l^, there is no universal 
expression for the conductance in the multichannel case; 
it depends on details of the electron distribution function 
in the leads. The latter may depend on specific geometry, 
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relaxation mechanisms, etc. BILP considered a par- 
ticular case of scattering-free leads connected to two ideal 
thermal reservoirs. Each reservoir is taken to be in ther- 
mal equilibrium with different chemical potentials on the 
left and right of the barrier. The distribution functions of 
the left and right moving carriers are strictly imposed by 
the reservoirs from which they originated. Thus on each 
side of the barrier left and right moving electrons occupy 
Fermi hemispheres of different radii. The potential drop 
across the barrier is taken to be the difference of effective 
chemical potentials of the leads, which are determined 
from the postulated condition of having equal number 
of particles above and holes below them. For the case 
of no channel mixing, Uj = USij and Vij — riSij, BILP 
obtained 

gBiLP = — r -I'l — rj— I (3) 

i 

where i runs over the open channels with velocity Vi. 
The same result was obtained earlier by Azbel us- 
ing a slightly different argumentation. For a single chan- 
nel Eq. (^) reduces to the old Landauer result, Eq. (||), 
whereas for a very high barrier, \ri \ k, 1, there is no differ- 
ence between the four-terminal result, Eq. (||), and the 
two-terminal one, Eq. (|^). 

In this paper we develop a logically clean formulation 
of the four-terminal conductance which avoids the prob- 
lematic introduction of different chemical potentials in 
various parts of the same system. Instead, we consider 
a pure Hamiltonian formulation of the problem with the 
external bias applied by a time dependent vector poten- 
tial through a macroscopically large loop containing a 
barrier. We solve self-consistent equations for the electric 
field on the Hartree (or RPA) level to determine the volt- 
age drop across the barrier. As a microscopic ingredient 
we need only an equilibrium linear conductance kernel, 
which is derived using the Kubo formula. Our formula- 
tion of the problem also allows us to tackle regimes not 
considered previously. As an example we consider a bar- 
rier embedded into a long multi-channel wire of uniform 
cross-section (a geometry different from that considered 
by BILP). If the wire length, L, exceeds vf/^jJ {vf is a 
Fermi velocity), the electron distribution function is a 
shifted (accelerated) Fermi sphere, rather than the BILP 
distribution. This results in a four-terminal conductance 
which is different from the BILP expression. For example 
in the low barrier limit, IrA <^ 1, we find 



(4) 



which is quite different from Eq. (^) in the same limit. 
For one thing, we do not expect discontinuities of con- 
ductance as function of the Fermi energy at the values 



where velocity of one of the channels is zero. In the oppo- 
site limit of a high barrier, \ti\ <C 1, we obtain Eq. (|l|), in 
agreement with BILP and other approaches. Another 
regime we are able to treat is a barrier embedded in 
a disordered wire with L ^ T', where is an elastic 
mean free path. In this case, the momentum equilibra- 
tion which takes place in the leads due to strong impurity 
scattering. We find that the same result, Eq. (^), holds 
for the case of disordered leads. That does not mean 
that the BILP expression, Eq. (^), is incorrect. We have 
recovered it, together with the conditions for its applica- 
bility, for a barrier embedded in a ballistic constriction, 
L < vf/(^,1'^\ between two wide reservoirs. Our proce- 
dure is quite different from that of BILP. In particular 
we do not require consideration of an unphysical, space 
dependent chemical potential — it may be defined using 
only standard equilibrium thermodynamics along with 
Poisson's and linear response equation. 

The outline of the paper is as follows: in Sec. || we for- 
mulate self-consistent equations for the induced electric 
field and charge density. In Sec. Ill we derive a micro- 



scopic expression for the non-local conductance in the 
case of the clean wire. We then solve the self-consistent 
field equations and obtain the four-terminal conductance 
in a number of regimes. Sec. |^ deals with a barrier with 
disordered leads. In Sec. ^ we treat an adiabatic con- 



Vl we discuss the 



striction geometry. Finally, in Sec. 
results, their domains of validity and possible generaliza- 
tions. 



II. SELF-CONSISTENT FIELD EQUATIONS 



Let us consider a multi-channel metallic wire of uni- 
form cross-section, 5, along the z-direction. We assume 
that in a vicinity of z = there is a potential barrier (or 
more generally some localized scattering region) across 
the wire. A weak uniform applied electric field, E'^^*{uj), 
with frequency lo creates a current density, which in the 
linear response regime has the form 



i{r,uj)= dr'g{r,r',uj) ■ E{r',uj) . 



(5) 



Here the integration runs over the volume of the wire 
and E{r',uj) is a total electric field at the point r' . The 
non-local linear conductivity kernel, g{r,r',uj), may be 
calculated using the Kubo formula. 

We assume that, due to the quasi-one-dimensional ge- 
ometry of the system, the electric field is practically inde- 
pendent of the transverse coordinates, E{r,u;) « E{z,uj) 
and only the total current I{z,uj) — dSj{r,uj) is mea- 
surable experimentally. The total current, I{z,uj), and 
the electric field, E{z,uj), are linearly connected through 
the conductance kernel, g(z,z',u!). Employing the Kubo 
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formula and disregarding the e-e interactions in 

the kernel g[z,z',(ju) (which is equivalent to the time- 
dependent Hartree approximation), one may express the 
non-local conductance as [13| 



/(z,z',w) 



(6) 



z, 2;')VzVz'G^(e_; z', z) , 



where AV^ B = AV^B - V^AB, e± = e ± u;/2, /(e±) 
is the Fermi function and G^'^^^ (e) is the retarded (ad- 
vanced) Green's function of the system at energy e in the 
transverse channel basis, 



G,,{e-z,z')^ dS dS'^*{x,y)G{e;r,r%{x',y'). 



Siz) S(z') 
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Here S,i{x,y) is the eigenfunction of the i-th transverse 
channel with energy e^. It will sometimes be convenient 
to divide the conductance into two parts. 



g{z, z', uj) = go{z - z' , uj) + gb{z, z' , uj) , 



(8) 



where go(z — z',lu) is the translationally invariant con- 
ductance of the wire without the barrier and gh{z, z' ,uj) 
is the contribution associated with the scattering by the 
barrier. 

The linear response relation, Eq. (H), together with the 
continuity and Poisson's equations, constitute a closed 
system of equations for three unknown quantities: the 
current I{z,uj), the electric field E^"'^{z,uj) and the 
charge density p{z, uj) all induced by the uniform external 
field, £;<=^*(w). 



dz'g{z,z',uj){E''''*{uj)+E'"'^{z',uj)), (9) 

(10) 
(11) 



/(z,c^) = 

dzl{z, uj) — iu!p{z, uj) , 
d,E''"^{z,uj) = 4ttp{z,uj)/S 



Eliminating I{z, uj) and p{z, uj), one obtains the following 
integral equation for E^"'^{z,uj) 



An 



d,E'^%z,uj) 
p^'Hz,uj), 



1 

iuj 



dz'dzg{z,z',uj)E'"'^{z',uj) 



(12) 



where p'--^") (z , uj) is the bare (unscreened) charge density 
produced by the external electric field near the barrier 



p(°)(z,uj)^^ fdz'd,g{z,z',uj)E'-\uj). (13) 



Since in the absence of the barrier no charge is induced 
in the leads, only the translationally non-invariant part 
of the conductance, gb{z, z' ,uj), contributes to p*^"-*. It 
is convenient to employ Fourier representation. Then, 
using Eq. (||), one may rewrite Eqs. (p^, (13) in the fol- 
lowing form 



'ujS 
in 
dq' 



-5o(<?,'^) 



E^-^{q,co) 



(14) 



2tt 

where e.g. 



g,{q, q', uj)E^"''{q', uj) = g,{q, 0, uj)E^^\uj) , 



gb{q, q', uj) = / dzdz' e'^'g^iz, z', c^)e-*«'^' . (15) 



Once Eq. (|lj) is solved, the induced voltage drop 
across the barrier is given by = J dzE^"'^{z,uj), where 
the limits of integration are taken to be much larger than 
any microscopic scale of the problem (see below) , but still 
much smaller than the length, C, of the ring which en- 
closes the driving time-dependent magnetic flux. In the 
Fourier representation one obtains 



V ^ E"''^{q~^0,ij). 



(16) 



where g ^ in such a way that <C (?o ^ uj/vp- The 
divergence- less part of the current, /q, originates from 
the convolution oi gQ{z — z') and E'^^'^: 



Io=go{q = 0,io)E--'{Lu) 



(17) 



All other terms in Eq. (^ describe currents localized near 
the barrier. As a result the four-terminal conductance, 
g = lo/V, may be expressed as 



g = go{0,uj) 



E'^'^^uj) 



E""^{q 0,w) 



(18) 



We shall proceed now by calculating the non-local con- 
ductance, g{z,z',uj), and then solving Eq. ( p^ for the 
various cases of practical interest. 



III. BARRIER WITH LONG CLEAN LEADS 

If there is no scattering of electrons in the leads, the 
Green's functions outside the barrier region may be ex- 
pressed as pj] 



G^Je; z, z') 



Uj{e)e'P^\'\+'P^\''\ , 



ipi\z\+ipj\z'\ 



zz' <0; 
zz' >0, 



(19) 



where pf/{2m) = ■mvf/2 = e — e^, etc. The transmission and reflection matrices of the barrier, tij{e) and ?'ij(e), obey 
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the unitarity condition 



EO^'<.(^)l' + K.WI'] =1- (20) 



Substituting Eq. ( p^ ) into Eq. (||), one can perform the energy integration assuming that <y(e) and rij{e) are slowly 
varying functions of energy on the scales of both temperature and frequency, w. As a result one obtains for the 
non-local conductance outside the barrier 



where — tijiep) and is the Fermi energy. Eq. (|2 
is the zeroth order term in the small parameter Lo/ep. 
Note, however, that even for small lo the combination 
Lozjvi is not necessarily small. 

Setting cj = in Eq. (Ell) leads to 



5 = g(z,z',c^ = 0) = — Trttt 

TT/l 



(22) 



The fact that at w = the conductance kernel is coor- 
dinate independent is a manifestation of the continuity 
equation. As a result the total current is divergence-less 
and given by Iq = g J dzE{z). The integral on the r.h.s. 
is a total voltage drop across the entire wire (and not 
just across the barrier). Therefore, the quantity g defined 
by Eq. ( p2[ ) represents the two-terminal conductance in 
agreement with the Landauer expression, Eq. (|l]). To 
analyze the four-terminal setup, one has to solve the 
self-consistent equation Eq. (|2|) and thus has to keep 
frequency the dependence of Eq. ( ^ ) . 

For simplicity we restrict ourself to the case of no chan- 
nel mixing inside the barrier region, i.e. tij — tiSij and 



~ ViSij, where, according to Eq. (Ed), |t 



1. 



In this case Eq. (^l|) may be simplified further (see Ref. 
i): 



(23) 



The two terms on the right hand side correspond to go 
and gb introduced in Eq. (||). In the Fourier representa- 
tion these terms take the form 



(26) 



iv,Gf{e;p)Gfie-p 



j tj - 1 + rj pp' tj-1- rj 



Here G'j^(e;p) = {e — ej — p^ / {2m) +ivi)^^ is the retarded 
Green's function of the j — th channel in the wire without 
a barrier (see Eq. (pj|)). 

By introducing the dimensionlcss momentum I = 
qvp/bj and normalized induced electric field E{1) = 
ujE""^{luj/vp,uj)/{vpE''''\uj)), Eq. (|l|) acquires the 
form 



E[l) 



dV 
2^ 



gb{l,l')E{l') ^ gb{l,Q) . 



(27) 



ci)^ EEE 8e'vp/{hS) = ujpT:hy{kj,S), 
Ane'^n/m being the plasma frequency. We have intro- 
duced the rescaled (frequency independent) conductance 



5o(0 



2v,. 
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2v, 



(28) 
(29) 



where Vi = Vi/vp. At small frequency, lo <^ ujp, one 
may neglect the first term on the l.h.s. of Eq. (p7|). The 
remaining equation is manifestly frequency independent. 
We shall, however, keep the frequency dependent term 
for the time being since it will allow us to discuss the 
characteristic length scales of the problem. We now pro- 
ceed to solve Eq. (E^ in some particular cases. 



ffo('7,'^) 



9b{q, q',uj) 



e2 ^ 



2v,; iuj 



El 



2viiw 2viiLu 



v^q^ — up' vfq'"^ — uj"^ 



(24) 



(25) 



These results may be derived directly in the momentum 
representation starting from the Kubo formula and em- 
ploying the following expression for the Green's functions 



A. Single— channel case 

First we apply the formalism to the well-known single- 
channel case. In this case the integral equation (j2^) has 
a separable kernel and may be easily solved. We neglect 
for a moment the frequency dependent term. Substitut- 
ing Eqs. (H) and (|2|) with vi = 1 into Eq. (|^), one 
obtains 



4 



dV E{1') 



= 1 . 



2i|r|2 J 27r [I'Y - 1 
The elementary solution of this equation is 



(30) 



(31) 



Employing Eq. (|18|), one immediately obtains the fa- 
mous Landauer expression, Eq. (|^) . This agrees with the 
conclusions obtained by Thouless ||] by imposing strict 
charge neutrality outside the barrier. Indeed, the fact 
that E^'"^{q) = const means that i;™'^(z) oc S{z). Thus, 
in this approximation there is no induced electric field 
and no induced charge density outside the barrier. To 
recover the shape of the charge distribution one must 
keep the frequency dependent term in Eq. (p7|). An ele- 
mentary calculation gives 



p(z, uj) cx sign(z) exp 




(32) 



showing that entire charge redistribution is confined to 
within the screening length, = vp/Cjp, near the bar- 
rier. Since we keep only the long wave-length compo- 
nents of the conductance kernel, Eq. (^) is valid only if 
the inequality k < kp is fulfilled (if not there are oscil- 
lations on the scale {2kp)~^). Although the initial set 
of equations, ([9|)-([Tl|), included the characteristic length 
vp/ui^ the expression (^) depends only on the much 
smaller scale, . This is a particular property of the 
single-channel case only. In the multi-channel case the 
length scale vp/u does not drop out of the solution of the 
self-consistent equations. This demands separate treat- 
ments for wires whose lengths are larger or smaller than 
Vp juj. 



If one keeps the frequency dependent term in Eq. ( 14 ) 
the resulting four-terminal conductance is 
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(33) 



At small frequency this expression describes the Lan- 
dauer resistance in parallel with an effective capacitance 



Co 



47r(2K-i) 



(34) 



which is the classical result for a plane capacitor of the 
area S and spacing 2k^^. In our approximations the 
screening length is assumed to be much larger than both 
the wavelength and the barrier width. If this is not the 
case, one has to keep the next order in Lo/ep as well as 
short wave-lengths in the expression for the conductance 
kernel, to obtain the effective capacitance. 



B. Equal velocity channels 

Another exactly solvable case, which will be used in 
Sec. ^ is that of N channels having the same velocity, 
Vi = 1; i = 1 . . . N . Calculations exactly parallel to the 
previous case lead to 



N 



(35) 



This result is in agreement with BILP, Eq. (|3|). If 
only M ^ N channels are open, ti = 0, |ri| ~ 1 for 
i = M -M . . . TV, Eq. (PI) simplifies to 



9 M 



(36) 



C. Weakly reflecting channels 

The case of a weakly reflecting barrier, |ri| <C 1 ; « = 
1 ... TV, may be considered by treating g;,(?, T) as a per- 
turbation in Eq. (p7|). Employing obvious operator no- 
tations and omitting the frequency dependent term, one 
may write the formal solution of Eq. (p^ as 



= {~%^9b + %^9b%^9b 



(37) 



where \5) = (1,0,0,...), where the first entry refers to 
I — component. Since is diagonal and may be eas- 
ily inverted, all terms in the expansion, Eq. (|37|), may be 
written in quadratures. We restrict ourself to the leading 
term only. In this approximation one obtains 



E{1) = i 



E4^|h|2(^2;2_i)-i 



(38) 



Substitution of Eq. (38) into Eq. ( [l8| ) leads to the result, 
announced in Sec. ^, 



1 2 



(39) 



Note that, apart from the single-channel case, the ap- 
proximate solution, Eq. (^), depends on the momentum, 
I, and hence on the z-coordinate. This means that strict 
charge neutrality does not occur in general. Rather, the 
resulting charge density exhibits a spatial modulation on 
the scale vp/uj. 
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D. Weakly transmitting channels 



We turn now to the high barrier case, ji^j <C 1 ; i = 
1 . . . A^. Employing |rip 
Eqs. (PI), m as 



= 1 — we rearrange 



where 

i ^ 



2v^ 



2i,. 



" 52/2 _i 52(;/)2_i 



(40) 

(41) 
(42) 



Here gi describes two disconnected half-wires, and gt - 
the perturbation due to the small tunneling transparency 
of the barrier. Since the conductance of the infinite bar- 
rier is zero, one expects that the operator gi has a zero 
eigenvalue. Indeed, it is easy to check that 



ffill) 



51^,0 = 0, 



(43) 



where |1) is the abstract vector in the momentum space 
given by (1,1,1...), whose entries refer to different val- 
ues of I. As a result, the operator gi is not invertible, 
which complicates the perturbation theory. To overcome 
this difficulty we pass to the basis of eigenstates of the 
operator gi , 



where a = 1,2,. 
takes the form 

{gt)ia'Ea' — - 



gi\a)^\a\a), (44) 
. and Ai = 0. In this basis Eq. (|27 



(.9fc)i ; 



Aa-Ba + {gt)aa'Ea' 



-{9b)a , 



a^l, 



(45) 
(46) 



where {gb)a — {(A9b\5), and summation over repeated in- 
dexes is assumed. Taking advantage of the smallness of 
gt |tj|2, one may solve these equations iteratively. Ne- 
glecting first gt on the l.h.s. of Eq. (|4^), one obtains in 
the zeroth approximation 



Ea 

El 



(gb)a 

(.9fe)i 

(at) 
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{9t)la' {gb)a' 



Xa' 



(47) 
(48) 



Repeatedly substituting the solution back to Eq. (46), 
one may obtain the higher order terms. In the leading 
(minus first) order in gt the solution has the form 



{gbh 

(5t)i: 



(49) 



To this order of accuracy on has to put jr^p = 1 in the 
numerator and therefore (56)1 = — (?o(0,w). By virtue of 
Eq. (^) (gt)ii = (.g)ii. After substitution into Eq. (|l|) 
one obtains for the four-terminal conductance 



9 = (5)1 



(50) 



This expression includes the exactly known |1) eigenfunc- 
tion only. A straightforward calculation of (5)11 gives the 
result, equivalent to the two-terminal expression. 



2 

e 



(51) 



This is the expected result, since for a high barrier, two 
and four-terminal measurements should yield the same 
result. Note, that to this order E{1) — const, which 
is consistent with strict charge neutrality. To evaluate 
corrections to the leading order, Eq. (pO|), one needs an 
explicit form of the eigenfunctions of the gi operator, 
Eq. (0), which is in general a hard problem. 



IV. BARRIER WITH DISORDERED LEADS 

We consider now the case, where the leads attached to 
the barrier contain weak elastic disorder. We restrict our- 
selves to uncorrelated short-range scatterers, with mean 
free path Z*^' and mean free time r related by l^^ — vpT. 
The leads are assumed to be much longer than the mean 
free path: L ^ I'^K As a result, the momentum distri- 
bution function far from the barrier is expected to be 
isotropic (spherical). If a current passes through the 
wire, the distribution function is a shifted Fermi-sphere. 
Our immediate goal is to calculate the disorder-averaged 
non-local conductance of such a system (leads with the 
barrier), substitute it to Eq. ( p^ ) and solve the latter. To 
this end one has to calculate the average product of two 
Green's functions, Eq. (^. The average single particle 
Green's function of the leads with the barrier is given by 
Eq. (^ ) where, in the disordered case 



Gf(e;p) 



1 



-p2/(2m) + i/(2T) 



(52) 



Employing the continuity relation, one may express the 
conductance kernel as 



9iq,q\^) 



qq' 



2TrSq,q' h iu}U{q, q' , uj) 



(53) 



The two terms in the square brackets represent the static 
and dynamic parts of the compressibility. In the diffusive 
approximation {kpl'^^ >> 1) the second term is given by 
the sequence of the bubble diagrams ||l^ , depicted in Fig. 
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|l|. Therefore the quantity Il{q,q',uj) may be found as a 
solution of the foUowing integral equation 



Uiq,q')=Tjiq,q') 



dq^ 



77(Q,g")n(g",g'^ 



(54) 



Here r^(7rwi) ^ is the inverse scattering amplitude and 

i 

r]{q,q',Lu) is the single bubble (see Fig. |l|) 



(55) 



where p± = p ± q/2, p'j. = p' ± q' /2 and the average 
Green's functions are given by Eqs. ( p6| ) and (^. Per- 
forming momentum summations, one obtains in the long 
wavelength limit {q,q' <ti kp, lo <^ ef) 

, v-^ T 1 — iUJT 



-y- 



nvi (1 — icur)^ + [viqrY 
\ri\'^Avfqq'T'^ 



[(1 - iujTf + {v,qTY][{l ~ iujtY + [v^q'Tf] ■ 

In the clean limit, lot ^ 1, the sequence is dominated 
by the single bubble resulting in 11 = 77. Under these 
conditions Eqs. (53) and ( p6| ) lead to the previous result 
given by Eqs. (|j)r(|5|). In the following we concentrate 
on the opposite, diffusive limit, vpqT <C 1, wt <^ 1. 




FIG. |l| Diagrammatic representation of the dynamic com- 
pressibility \l{q,q' denoted by the shaded bubble; the 
quantity r){q,q' ,uj) is denoted by the empty bubble. The full 
lines represent electron Green's function, Eqs. (]2^), ( js^ ) and 
the full dot represents a point-like elastic scatterer with the 
amplitude {tv)~^ . 

In the absence of the barrier, — 0, and Eq. ( |5^ ) may 
be easily solved, resulting in 



Ii{q,q',uj) 



Dq^ — ioj 



Iloiq,Lu)2TTS,^g, , (57) 



where the single particle density of states (per spin), v, 

de: 



and the diffusion constant, D, are defined as 



D 



(58) 



Substituting Eq. ( p7D into Eq. (|5^) , one obtains the stan- 
dard diffusive expression for the conductance of a disor- 
dered wire 



90 



2vAujT 



nh Dq^ — iuj 



(59) 



Next we shall look for gb{q,q',oj) and then for solutions 
of the self-consistent equation (|l^) in some particular 
cases. 



A. Equal velocity channels 

In the case of N channels with the same velocity, 
Vi — V, i — 1 . . . N, the integral equation ( p^ has a sepa- 
rable kernel and thus may be solved exactly. After simple 
algebra one obtains for the barrier-induced part of the 
conductance 



gb{q,q',uj) = 



2viLUT 2viujT 



Dq^ — iuj Dq'- 



(60) 



This expression is valid for ut <C 1; however, we have 
retained the frequency dependent term in the denomina- 
tor since the condition ^/uPF ^ l^jP/kjP is not assumed. 
Substituting Eqs. (|5|), ( |60| ) into the self-consistent equa- 
tion for the induced electric field, Eq. (pH), and solving 
the latter, one obtains, for a; <C cjp. 



E'^'^iq) = 2vT 



(61) 



Finally, employing Eqs. (|lj) and (^9|), yields the four- 
terminal conductance 



N 

j 



(62) 



El 



in agreement with a clean case, cf. Eq. (pS]). Although 
for the sake of compactness we wrote all intermediate 
expressions for the small frequency limit, lot <C 1, one 
may actually perform all the calculations and obtain the 
final result, Eq. (|6^), for any lot (the only limitation is 
CO <^ Cop). 
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B. Weakly reflecting channels 

One may formally solve the integral equation ( |5^ ) em- 
ploying perturbation theory with the small parameter, 
|rip <C 1. The first order correction to Il{q,q',uj) for 
wr <C 1 is 



SU{q,q',u;) 



T 

2^ 



El 



2viq 



2v,q' 



Dq^ — iuj Dq'- 



(63) 



Employing Eq. (p3|), one obtains in the same order 
9b{q,q',uj) = - 



2 -. (64) 



T:h ^ D(P- — iu Dq'"^ — iu. 

i 

Note that Eqs. (^), ( |64| ) are almost exact analogs of 
Eqs. (p^), (^), where the ballistic propagators are re- 
placed by diffusive ones. The important difference, how- 
ever, is that the validity of Eqs. (p^), ( ^ ) is not restricted 
to small reflection coefficients, r^. On the other hand, in 
the diffusive problem an electron can bounce between the 
barrier and impurities, and thus the exact conductance 
contains higher powers of Ir^p. Eqs. ( ^9[ ) and ( |6^ are 
only the first two terms in the infinite series. 

We now substitute Eqs. (||), (||) into Eq. ( P) and 
solve the integral equation. Since the kernel, Eq. (^4|), is 
separable, its solution is elementary. The result for the 
induced electric field in the zero frequency limit is 



(65) 



Employing Eq. ( |18| ) and ([59D, one obtains the same re- 
sult for the four-terminals conductance, g, as in the clean 
case - Eq. (Q). This statement is actually valid for any 
value of uiT, provided |rip <C 1 and a; <C Wp. 



C. Weakly transmitting channels 



For a completely reflecting barrier, jr^ 
easily check the following equality 



1, one may 



2^ 



2t:6, 



1 



9:9 



VT 



7^" 



With its help and employing Eqs. (53), (& 



with 



1, one may prove that 



2tt 



9'"' ^{q,q',^) 







(66) 
and ^) 

(67) 



for any q and uj. This identity simply reflects the fact 
that an electric field localized under a completely reflect- 
ing barrier cannot induce any current. Thus, to solve the 
self-consistent field equation for \ti\ <C 1 one has to invert 



an operator having one almost vanishing eigenvalue. Fol- 
lowing the method described in section III D| , one obtains 
for the the four-terminal conductance 



(5) 
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(68) 



where {g)ii = (27r)~^/ dqdq'g{q,q'). Calculation of this 
matrix element is not as straightforward as in the clean 
case. To compute it one has first to solve the integral 
equation (|5j). Simple algebra reduces that problem to 
the inversion of the operator which differs by a small fac- 
tor (~ from the one written in the square brackets 
on the left hand side of Eq. (|66|). According to Eq. ( |6^ ) 
the latter has exactly zero eigenvalue, with the corre- 
sponding eigenvector cx 1/q'. One has to employ once 
again the method described in section [II D| to invert an 
operator with one small eigenvalue. As a result of this 
procedure one obtains Eq. (EV) , the same as in clean case. 



V. BARRIER IN A BALLISTIC CONSTRICTION 



We now consider a barrier embedded into a ballistic 
adiabatic constriction between the two wide reservoirs 
(see Fig. ^a). We shall assume that the length of the 
constriction, L, satisfies the conditions L < vp/oj and 
L < l^K This is the geometry considered by BILP |Tc| ] 
and others (see Ref. |^ and references therein). Since the 
cross-sectional area of the constriction, S{z), is a smooth 
function of z, one may write the wave functions in the 
adiabatic approximation p^ ] 



(69) 



where ^i{x,y; S{z)) is a transverse wave function of the 
z-th channel with an eigenenergy ei(S{z)) = ei{z). The 
longitudinal function, ipi{z), satisfies a one dimensional 
Shrodinger equation 

[-(2m)-ia2 + g^(^) + Vbiz)]Mz) = ^M^) , (70) 

where Vfc(z) is the localized tunneling barrier potential. 
Thus the problem is reduced to a one dimensional effec- 
tive tunneling problem described by Eq. (|7^). (See also 
Fig. H) If the constriction is narrow enough it closes 
most of the channels, since in these channels ei{z) > ep 
for small enough z. Hereafter we shall assume that this 
is the case. Even in the open channels there is a partial 
reflection of electrons due to the presence of the barrier, 
which is characterized by the barrier's reflection coeffi- 
cients, Ti. For z and z' outside the barrier region, by 
employing Eq. ^ and expressing the Green's functions 
in the WKB approximation, one flnds for the conduc- 
tance kernel (cf. with Eq. (p^) ) 
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z 






f dz 


exp < 


lUJ 


J v,{z) 










z' 



(71) 



\n\ exp < lu 



dz 

Vt{z) 



+ iuj 



dz 

Vi{z) 




• A .C 



• D .B 



where 



Vi{z) 



2{eF-e^iz)) 



(72) 



and the summation runs over the open channels only. 
This expression will be used to find the four-terminal 
conductance of the barrier, measured by probing the volt- 
age drop across it. (In very narrow energy intervals, when 
one of the channels is extremely near its transmission 
threshold it may not be tractable by WKB.) 

Let us first imagine solving the problem on the large 
length scale, |z| ^ L where the entire constriction 
may be treated as a single localized scatterer. The re- 
sults of previous sections are directly applicable to such 



a problem (see e.g. Sec. IIIB). One finds that at suffi- 



ciently low frequencies the induced electric field, E , is 
confined to the constriction region. A line integral over 
this field is a total voltage drop across the constriction 



(b) 




E'"'(z) 




(d) 



dzE"'"^{z,oj) = Vo 



(73) 



where the integral effectively runs over the region be- 
tween the two points A and B (see Fig. ^) located 
outside the constriction at a distance which is large com- 
pared with a bulk screening length (but still much smaller 
than L). Since only a few channels are open, the current 
is determined by Eq. (|36|) and given by 



Vn 



(74) 



There are corrections to this expression of the order of 
a ratio of number of open channels to a total number of 
channels Only the open channels contribute to the 
sum on the right hand side; the transmission coefficients 
of the open channels are those of the localized barrier. 



FIG. g (a) Schematic drawing of the adiabatic constriction 
with the barrier, (b) The induced electric field, E^"'^{z), given 
by Eqs. (|76|), ([ts]). (c) The scalar potential which corresponds 
to the induced electric field, (z). (d) The efi'ective poten- 
tial ei[z) + Vb[z) for the open (full line) and closed (dashed 
line) channels. The dotted line represents the Fermi energy. 



To determine the voltage drop across the barrier it- 
self one must determine the structure of E'^^'^{z,uj) on a 
small length scale, which turns out to be the screening 
length, K^^(<C L). The induced field E'^"'^ is comprised 
by the two components, i?™'^ = iJ™'^ -|- iJ™'^, created by 
charge densities induced at spatially well separated lo- 
cations. El""^ is created by charges accumulated in the 
regions of the narrowing of the constriction. This field 
component is practically constant between points C an D 
in the constriction (see Fig. ^ja). Its contribution to the 
voltage drop measured across the barrier is small in the 
parameter Zy/L <C 1, where 2zy is a distance between 
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the voltage probes. The other component, £^2"'^, is in- 
duced by charges accumulated near the localized barrier. 
It is sharply peaked in the barrier region and is primarily 
responsible for the voltage measured across the barrier. 
The induced field, w), inside the constriction re- 

gion is given by a solution of the self-consistent equation 
(p^), with the conductance kernel Eq. ([7l|). For low fre- 
quencies, (ujL/vf) ^ 1, one may expand the exponents 
in Eq. (^TJ) it terms of this small parameter and obtain 
to leading order 



hsjz) 



E 



sign(2:)Vo : 



(75) 



where Eq. ( |73| ) was employed on the right hand side. In 
the region where the effective potential is "fiat" (see Fig. 
||d) and where Vi(z) k, Vi = const and S{z) k, S ~ const, 
the solution of Eq. (|75|) is 



2,,-i 



£;"'^(z) = const + K, -p 

i 

where the screening length is 



exp{— |z|k} 



(76) 



VF 



(77) 



The constant on the right hand side of Eq. (|7^) represents 
El"''-, the component created by the distant charges of 
the constriction's narrowing. Its value may be estimated 
from Eqs. ^ and {M) as 



Vn 



Er « const « -f 



(78) 



The voltage drop across the barrier, measured by the 
probes located just outside the barrier at points E and 
F (see Fig. ^) with coordinates z = ±Zy, is given by 
V = dzE'"''. Employing Eq. (|7^), one finds 



V = Vo 



i 



(79) 



Therefore in the asymptotic regime k ^ ^ z„ <C i, one 
finds for the four-terminal conductance g = Iq/V, 



^ i i 



(80) 



in agreement with BILP |10|. We therefore confirm 
their result for the barrier embedded into the constric- 
tion attached to the wide reservoirs if the inequalities 
<C i < Vi/u), l"^' are satisfied. 
Eq. ( |80| ) predicts a sudden drop of the conductance if 
one (or several) 0. The physical reason for such 

behavior is an enhanced screening by the channel(s) with 
the small velocity and hence a large density of states. We 
stress that although a non-monotonic dependence of the 
conductance on the chemical potential may indeed take 
place, there is no actual discontinuity when a channel 
opens. The reason is that the above calculations result- 
ing m Eq. (1^ lose their validity if < h/{mL) ^ (or 
Win < Lbj). The two key approximations simultaneously 
go wrong in this case: (i) the long-wavelength limit for 
the non-local conductance, which neglects Friedel oscilla- 
tions at q — 2kF and (ii) the WKB approximation, which 
neglects reflection by £10(^)1 the effective barrier due to 
the constriction. If these two effects are properly taken 
into account the continuity of conductance is restored, 
however its value for < < Ti/ {mL) depends both on 
the position of the probes and the specific shape of the 
constriction. 

Also, at the beginning of this Section, we assumed that 
if a channel is "closed", i.e. Vi ivi, the corresponding 
current vanishes. In fact this requires negligible tunnel- 
ing or \via \ > h/{mL). Thus Eq. ( |80| ) also fails just below 
an energy threshold. 

Finally we address the frequency dependence of the 
four-terminal conductance. To this end we expand 
Eq. ( [7l| ) to the next order in LoL/vp and substitute into 
the self-consistent equation (p^). This gives the fre- 
quency dependent correction to the induced electric field 



7ind(. . I exp{~|z|K} 



2k- 



dz\z\E''"'{z) , 



(81) 



where £^'"'^(2) is the zero frequency field given by 
Eqs. (|^), (II). According to Eqs. (|^), (|3|) the rela- 
tion between the divergence-less current and the total 
voltage drop, Eq. ([z^, should be modified to include 
the classical capacitance of the constriction, Cq. The re- 
sulting four-terminal conductance acquires a capacity- 
like frequency dependent correction of the following form 
5g{uj) = —iuj{Co + C), where 



Jdz\z\E''"''^{z) 



Zlnl^v-' !dzE^"d{z) 



(82) 



The precise value of the last factor in this expression 
depends on E"'''^{z) on the scale z ^ L and hence on 
the shape of the constriction. It may be estimated using 
Eqs. (|7|). Provided that {Lk)~'^ < \U\'^, which is 
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valid for all barriers except those with almost complete 
reflection, one finds for the additional capacitance 



^ = E l^^l' '-^^ ^ (83) 



L \rt\'^v. 



2 -1 



0) 



where 7 = 0(1). In the opposite limit |tip ^ (-^'«) ^ 
one obtains 



(84) 



2„-l 



III 



II 



L 



We stress the peculiar dependence of the additional ca- 
pacitance on channel velocities and in particular a sharp 
rise at the point of a new channel opening. (As was 
mentioned above, there is no actual discontinuity, since 
Eqs. (|83|), ( p4| ) are not applicable once < fi/{mL).) 



VI. DISCUSSION OF THE RESULTS 



We have presented an approach to the Landauer four- 
terminal tunneling conductance, which does not assume 
the presence of two reservoirs which maintain different 
chemical potentials. Our theory is based on the solu- 
tion of a self-consistent set of equations involving the 
non-local conductance, 5(2, z' , uj), the tunneling current, 
induced charge density and electric field. 

There are three different regimes (see Fig. ^), defined 
by the length of the leads, i, the frequency, w, and the 
rate of elastic scattering, r^^ (or the elastic mean free 
path, l"^^ — vft). The first "high" frequency regime 
denoted by I in Fig. ^ is defined by the conditions 
uj > maxjr^^, tii?/L}. We have treated it in section 
The second denoted by II is that of "disordered" leads. 
It is defined by the conditions L > F'; < and was 
treated in section Finally the third one denoted by 
III and defined hy L < Z*^'; w < vp/L is the regime of 
"ballistic" motion. It was considered in section M. 



FIG. ^ Various physical regimes; area I corresponds to a 
"high" frequency regime, II to a "disordered leads" and III to 
a "ballistic constriction" geometry. Here L is the length of the 
leads attached to the barrier, u - the frequency of the driving 
field, and are elastic mean free path and scattering 
rate correspondingly. 

A few remarks about the relationship of our work in 
the ballistic regime with a barrier, in the geometry of Fig. 
Ija, and the work of BILP |jl^: First, the final result (our 
Eq. (|8^) ) is the same in both treatments. However: in 
BILP the occupation of electronic states in the presence 
of a current is different in the regions to the left and to 
the right of the barrier. In our work the presence of a 
current does not modify the occupation of the electronic 
states (which extend from the left through the barrier to 
the right), but does modify their wave-functions - in ac- 
cordance with standard Kubo transport theory. BILP's 
chemical potential differences ^1 — ^2 and ^a — IJ^b are re- 
spectively equal to the electrostatic potential differences 
Vb and V of Fig. We have not been able to derive 
the very simple picture of BILP from our considerations. 

In the regimes I and II we have obtained results which 
are qualitatively different from ballistic regime III. The 
most striking difference occurs for the weakly reflecting 
barrier (cf. Eqs. (^ and (^)). In particular, for long 
and/or disordered leads the conductance does not exhibit 
discontinuities when the velocity of one of the channels 
vanishes. The physical reason for these differences lies in 
the difference of momentum distribution functions in the 
diffusive and high frequency regimes compared to that 
in the ballistic regime. In the opposite limit of weak 
transmission one finds, to leading order, the same result 
for all three regimes, which is identical with Landauer's 
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two-terminal expression, Eq. (|l|). We believe that the 
predicted qualitatively different dependence on the chan- 
nel velocities (and hence on the chemical potential) for 
diffusive and ballistic leads may be checked experimen- 
tally. 

We also predict a peculiar dependence of the capaci- 
tance on the chemical potential, Eqs. (|8^), (|^, which, 
we hope, may also be checked experimentally. 
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